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Abstract Employing an approach based on the Green functions of Landau-gauge QCD, some selected results 
from a calculation of meson and baryon properties are presented. A rainbow-ladder truncation to the quark 
Dyson- Schwinger equation is used to arrive at a unified description of mesons and baryons by solving Bethe- 
Salpeter and covariant Faddeev equations, respectively. 
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1 Motivation: 

Why Functional Approaches to 
QCD? 

The central aim of the studies reported here is 
to develop a QCD-based description of the structure 
of hadrons in terms of quarks and gluons. Theoreti- 
cal issues such as confinement, dynamical breaking 
of chiral symmetry and the formation of relativis- 
tic bound states can be understood and related to 
the properties of QCD's Green functions. Significant 
progress within functional approaches has recently 
been achieved especially in Landau gauge, see e.g. 
refs. ti'— 'iJ. Hereby it should be noted that differ- 
ent pictures of confinement are not necessarily mu- 
tually exclusive but may describe different facets of 
this phenomonen, see ^ and references therein. The 
advantage of employing Green function methods is 
given by the fact that these Green functions provide 
the input for QCD's bound state equations which, in 
turn, can be used to calculate hadron properties. The 
Green functions of elementary fields are determined 
by Dyson-Schwinger equations, Exact Renormaliza- 
tion Group equations, and/or lattice calculations. It 
should be noted that dynamical breaking of chiral 
symmetry leads to the generation of quark masses 
and scalar quark-gluon interactions 

Any numerical solution of QCD's bound state 
equations, i.e. the Bethe-Salpeter equation for 
mesons and the Faddeev equation for baryons, re- 
quires Green functions of quarks and gluons as input, 



which in turn necessitates a truncation of the cor- 
responding functional equations. The results which 
will be presented in this status report build upon a 
rainbow-ladder truncation, i.e. a dressed-gluon ex- 
change, for mesons and baryons in a unified ap- 
proach ILAAifl'iiJ. A symmetry-preserving exten- 
sion to more realistic kernels is on its way for mesons, 
see anc l references therein. For baryons as 
three-quark states, the covariant Faddeev equation 
in rainbow-ladder truncation has been solved only re- 
cently 112. 

2 Bound state equations and 
rainbow-ladder truncation 

Hadrons, being bound states of quarks, appear as 
poles in the qq and qqq scattering matrices. Hadron 
properties can be extracted upon solving bound-state 
equations which are valid at these poles and need the 
elementary QCD Green functions as input: 



(1) 



Here ^ is the bound-state amplitude defined on the 



hadron's mass shell, and K 



C-2) 



K (2) SS, K, 



(3) 



S S S are the renormalization-group invariant 
products of the respective kernel Ki n \ with two or 
three dressed quark propagators. 

These fully Poincare-covariant bound-state equa- 
tions provide a tool to calculate experimentally acces- 
sible hadron observables, e.g. meson and baryon mass 
spectra, decay constants, scattering processes, and 
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Fig. 1. The quark-diquark BSE in pictorial form. 



electromagnetic properties such as form factors, mag- 
netic moments and charge radii (see e.g. Refs. l&^J 
and references therein). To proceed with the nu- 
merical solution of a covariant bound-state equation 
like the Bethe-Salpeter or the Faddeev equation, one 
needs to specify all ingredients: the interaction ker- 
nels and the dressed quark propagator which ap- 
pear in these integral equations. These Green func- 
tions can be determined either from functional equa- 
tions and/or lattice calculations. In addition, they 
are related via Slavnov- Taylor and Ward-Takahashi 
identities. With respect to the pion, the would-be 
Goldston boson of chiral symmetry, the axial-vector 
Ward-Takahashi identity plays a special role. If it is 
satisfied by the interaction kernels in related equa- 
tions, this correctly implements chiral symmetry and 
its dynamical breaking, leading e.g. to a general- 
ized Gell-Mann-Oakes-Renner relation which is valid 
for all pseudoscalar mesons and all current-quark 
masses liS-JS. In particular, the pion is massless in 
the chiral limit independent of the details of the in- 
teraction. 

Although the structure of the kernel is restricted 
by the above mentioned identities, to construct a 
symmetry-preserving truncation of the bound-state 
equations is a non-trivial task. The rainbow-ladder 
truncation provides such a symmetry-preserving 
scheme. Instead of interaction kernels being func- 
tions of the quark and gluon momenta one needs 
a functional form for the interaction as a function 
of the gluon momentum only. The non-perturbative 
dressing of the gluon propagator and the quark-gluon 
vertex are "absorbed" into an effective coupling for 
which we adopt a widely-employed ansatz liiiS. It 
reproduces the logarithmic decrease of QCD's one- 
loop perturbative running coupling at high momenta, 
and its infrared contribution is parametrized by an 
infrared scale and a dimensionless parameter (for de- 
tails, see ^). Especially, it yields the non-pertur- 
bative enhancement at intermediate gluon momenta 
necessary to generate dynamical chiral symmetry 
breaking, and hence a constituent-quark mass scale. 



We emphasize here that the input of the Faddeev 
equation is completely specified at the beginning with 
all parameters already fixed from a few meson prop- 
erties. In addition, we want to note that the resulting 
quark propagator is in very good agreement with cor- 
responding results of fully coupled Dyson-Schwinger, 
resp., Functional Renormalization Group equations 
and of lattice calculations. 

3 Baryons: Quark-diquark model 

As the amplitudes resulting from a fully relativis- 
tic Faddeev equation are of high complexity it is natu- 
ral to look first for an approximation keeping Poincare 
invariance and the relation to QCD Green functions. 
A well-studied approximation of this type is based 
on using diquark correlations as effective degrees of 
freedom. 

The motivation for studying diquark correlations 
has been the observation of a strong attraction in 
the SU(3)c antitriplet qq channel, e.g. in lattice 
[19,20, ^m^m, and Bethe-Salpeter studies. 
Such an attraction has also been proposed to explain 
missing exotic states in the hadron spectrum and the 
masses of light scalar mesons ISL^SJ. Further support 
for the diquark concept has been provided by a study 
of diquark confinement in Coulomb-gauge QCD [29 ' . 

These arguments lead to the assumption that cor- 
relations between two quarks provide the dominant 
attraction not only in meson but also in baryon chan- 
nels. Consequently, the quark-diquark model traces 
the nucleon's binding to the intrinsic formation of col- 
ored scalar- and axialvector diquarks. It treats such 
two-quark correlations as a separable pole sum in the 
qq scattering matrix which leads to a description of 
baryons as bound states of effective quarks and di- 
quarks. In this way the covariant Faddeev equation 
(see Section [4]) is simplified to a two-body bound- 
state equation (see Fig. [IJ while full Poincare covari- 
ance is maintained. 

Nucleon and A properties have been studied in a 
quark-diquark model with parametrized ingredients, 
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Fig. 2. Faddeev equation J2J in rainbow-ladder truncation. 



see e.g. \2&Mi£2i£2i an( j references therein. The ap- 
proach was subsequently extended to determine the 
dynamics of + and 1 + diquarks from their underly- 
ing quark and gluon constituents 1^ — ' — ' — ' ^ where 
parametrizations for the diquark amplitudes were re- 
moved and replaced by solutions of the corresponding 
diquark Bethe-Salpeter equations. 

The identification of colored diquarks as poles 
in the qq scattering matrix is possible within a 
rainbow-ladder truncation: one obtains timelike + , 
1 + , . . . diquark poles. The corresponding mass scales 
play an important role in the description of light 
baryons. These poles, however, correspond to un- 
physical asymptotic states and disappear from the 
spectrum when going beyond rainbow-ladder. The 
obtained rainbow-ladder diquark masses exhibit large 
sensitivities to an effective width parameter of the 
rainbow-ladder interaction kernel, a feature which has 
previously been observed in Ref . ^ . 

The results for nucleon and A masses employ- 
ing the quark-diquark bound state equation (Fig. [1]) 
are shown in the left panel of Fig. [3] together with 
a selection of lattice results. The corresponding ab- 
scissa values mjj., as well as m p , are obtained from the 
pseudoscalar and vector-meson Bethe-Salpeter equa- 
tions in one consistent calculation. The findings are 
qualitatively similar to those for m p : setup A, where 
the coupling strength is adjusted to the experimental 
value of f w , agrees with the lattice data. This be- 
haviour can be understood in light of a recent study of 
corrections beyond rainbow-ladder truncation which 
suggests a near cancellation in the p-meson of pionic 
effects and non-resonant corrections from the quark- 
gluon vertex Setup B provides a description of 
a quark core which on purpose overestimates the ex- 
perimental values (see the discussion below) while it 
approaches the lattice results at larger quark masses. 

4 Baryons: 

covariant Faddeev equation 

A formalism to treat the three-body bound-state 
problem through the analogue of the two-body Bcthc- 



Salpeter equation was formulated in Refs. 125? see 
122 for an overview. The corresponding covariant Fad- 
deev equation describes the baryon as a bound state 
of three spin- 1 /^ valence quarks (cf. Fig. O: 

3 

where the interaction kernel comprises a three- 
quark irreducible contribution and the sum of per- 
muted two-quark kernels whose quark-antiquark ana- 
logues appear in a meson Bethe-Salpeter equation. 
The subscript a denotes the respective accompanying 
spectator quark. 

The QCD Green functions once more provide a 
way to embed this equation in a consistent quantum- 
field theoretical setup. Its dynamical ingredients can 
then be treated in perfect correspondence with stud- 
ies of quark and meson properties. A solution of 
the equation relies upon knowledge of the dressed 
quark propagator and the three-quark kernel. A nu- 
merically solvable form of this equation is based on 
the specification and decomposition of the Poincare- 
covariant baryon amplitude. 

Recently a fully Poincare-covariant computation 
of the nucleon's Faddeev amplitude was performed 
in such a unified treatment U&^SJ. The assumption 
of dominant quark-quark correlations motivates the 
omission of the three-body irreducible contribution 
from the full three-quark kernel in Eq. ([2]). The 
numerical solution of the Faddeev equation is then 
performed after truncating the interaction kernel to 
a dressed gluon-ladder exchange between any two 
quarks, thereby making a direct comparison with cor- 
responding meson studies meaningful. Further in- 
sight into the structure of baryons is provided by 
contrasting the results to those of investigations of 
baryons in the quark-diquark model. 

The resulting nucleon mass M N and its evolution 
with is plotted in the right panel of Fig. [3] A com- 
parison to the consistently obtained quark-diquark 
model result exhibits a discrepancy of only ~ 5%. 
This surprising and reassuring result indicates that a 
description of the nucleon as a superposition of scalar 
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Fig. 3. (right panel adapted from Ref.^) Evolution with m% of m p , Mjv and Ma compared to lattice 
data The left panel shows the quark-diquark model results for nucleon and A masses ^'^J; the right panel 
compares the quark-diquark and Faddeev results for the nucleon mass. Dashed and dashed-dotted lines 
correspond to setup A (i.e. f n is fixed as input); the solid line for m p and the bands for Mjv and Ma are the 
results of setup B (advertent inflation of hadronic observables to leave room for pionic corrections) , where the 
variation w.r.t. to the interaction width parameter is explicitly taken into account to provide error estimates. 
Dots denote the experimental values. 



and axial-vector diquark correlations that interact 
with the remaining quark provides a close approxima- 
tion to the consistent three-quark nucleon amplitude. 

By construction, the rainbow-ladder truncation 
omits all mesonic dressings of baryons. As pionic ef- 
fects are certainly important for all hadrons this leads 
to the idea of changing the effective interaction such 
to advertently inflate the p mass (Set B in Fig. [3]) and 
thus to leave room for pionic corrections at small pion 
masses. It can be seen that such a treatment effects 
also the baryon masses in a way which provides room 
for further attraction by pions. 

5 Nucleon Electromagnetic 
Form Factors 

In order to calculate the nucleon's electromag- 
netic form factors in the given approach one has to 
relate the nucleon's electromagnetic current to the 
underlying description of the nucleon as a compos- 
ite object. To this end the baryon must be resolved 
into its constituents to each of which the current can 
couple. A systematic procedure for the construction 
of a hadron-photon vertex based on electromagnetic 
gauge invariance 121'— '^2J yields for the current oper- 
ator 

,r = -% (Gff + Go K» Go) tf< , (3) 



where G denotes the product of dressed propagators 
and K the kernel which appears in the respective 
bound-state equation. In the quark-diquark model, 
the incoming and outgoing baryon states are de- 
scribed by quark-diquark amplitudes <!>;, <&/. Upon 
interaction with the external current, the baryon is 
resolved into its constituents: quark and diquark and 
the interaction between them. To each of these ele- 
ments the current couples t2i^J. 

Previous nucleon form factor studies performed 
within the quark-diquark model ^2? £L — -Li- i£? -21 
share some common caveats. First, as already in- 
dicated above pionic contributions play an important 
role in the low-energy and small-quark mass behav- 
ior of the nucleon's electromagnetic structure. Such 
effects are not included in a quark-diquark 'core' and 
must be added on top of it 122? i2? i2? £SJ . Second, access 
to the large-Q 2 region and thereby to the truly per- 
turbative domain is so far only feasible upon imple- 
menting pole-free model propagators which, in turn, 
exhibit essential singularities at timelike infinity, cf. 
ref. . The problem of taking into account the real- 
istic analytic structure of the quark propagator is not 
of fundamental concern; it merely awaits a thorough 
numerical treatment. Third, the quark-mass depen- 
dence of magnetic moments and charge radii, while 
emerging naturally in lattice calculations, is practi- 
cally inaccessible in a quark-diquark model due to 



No. X 



R. Alkofer et al: Hadron properties from QCD bound-state equations 



5 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 





(r/) 2 [fm 2 ] 















1.0 



0.6 



0.4 



0.2 



0.0 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 

m n [GeV] m n [GeV] 

Fig. 4. (adapted from Ref.^) Quark-diquark model results for the squared isovector radii (corresponding to 
the Dirac and Pauli form factors Fi j2 = Ff 2 — F"^) in setup B, compared to lattice results^'— ^J. 



the unknown mass dependence of the modeled ingre- 
dients. 

The rainbow-ladder based quark-diquark ap- 
proach described in Section [3] removes the latter ob- 
stacle. Upon resolving the diquarks' substructure, 
the form factors are immediately related to the pa- 
rameters in the effective quark-gluon coupling a(k 2 ), 
in particular: its quark-mass dependent coupling 
strength. Here we restrict ourselves to the quark 
'core model' which represents a quark-diquark core 
that needs to be dressed by meson-cloud effects. A 
comparison of the core's static properties with lattice 
results is appropriate at larger quark masses; form 
factors depending on the photon momentum may be 
compared to experiment at Q 2 > 2 GeV 2 where pion- 
cloud effects are diminished. 

The nucleon's Dirac and Pauli radii r x and r 2 are 
shown in Fig.H] and follow a similar pattern as the 
pion charge radius they are weakly dependent 
on the width parameter of the interaction and agree 
with lattice data at larger quark masses where the 
'quark core' becomes the baryon. A natural feature 
of a quark-diquark model is the negativity of F"(Q 2 ). 
The presence of an axial-vector dd diquark correla- 
tion centers the d-quark in the neutron and induces 



r d [ill 
1 



Only the components transverse to the photon 
momentum, i.e. those not constrained by current con- 
servation, determine the physical form factor content. 
The respective contribution to the quark-photon ver- 
tex is known from its inhomogeneous Bethe-Salpeter 
solution 1£2J an d includes a p-meson pole which 
amounts to ~ 50% of both pion and nucleon squared 
charge radii. The available information on those parts 



at larger Q 2 , especially for the diquark-photon in- 
gredients, is limited within the scope of the current 
approach; however they are mandatory to enable a 
realistic Q 2 -evolution of the form factors and the pro- 
ton's form factor ratio n v G p E /G p M It is notewor- 
thy that such uncertainties will be removed upon ex- 
tending the the quark-diquark approach to a three- 
body framework: in a rainbow-ladder truncation, the 
analogue of Eq. ((3J only involves the dressed-quark 
photon vertex and a dressed gluon propagator, i.e., 
quantities which are known and need not be modeled. 

6 Summary and Outlook 

In Landau gauge QCD, the Green functions of 
quarks and gluons are sufficiently well determined 
such that they can be used as a reliable input in rela- 
tivistic bound state equations. These provide then 
a unified approach to mesons and baryons within 
a quantum-field-theoretical framework. Within the 
rainbow-ladder truncation meson observables and the 
nucleon mass have been calculated. These recent re- 
sults mark a milestone in the effort to determine nu- 
cleon observables in a 'first-principle' functional ap- 
proach to continuum QCD. 

Electromagnetic properties of the nucleon and 
other baryons are currently investigated within this 
truncation scheme. The extensions of recent stud- 
ies of meson properties beyond rainbow-ladder ™ to 
the covariant Faddeev equations will improve on the 
now existing unified approach to meson and baryon 
properties. 
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